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Born-lnfeld non-linear electrodynamics was introduced to render the self energy of a point par- 
ticle finite. It has recently been revived as a field theory for branes and strings. We quantize 
this theory on a Euclidean space-time lattice, using Metropolis Monte-Carlo simulations to mea- 
sure the properties of the quantum field theory. Liischer-Weisz methods are used to measure the 
electromagnetic fields from a static point charge. The D field from a point charge appears to 
be identical to that for the normal Maxwell Lagrangian. The E field is enhanced by quantum 
fluctuations, and shows short distance screening as it does in the classical theory. 



XXIIIrd International Symposium on Lattice Field Theory 

25-30 July 2005 

Trinity College, Dublin, Ireland 



* Speaker. 

'This work was supported by the U.S. Department of Energy, Division of High Energy Physics, Contract 
W-31-109-ENG-38. 

* Supported in part by NSF grant NSF PHY03-04252. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it/ 



Lattice simulations of Born-Infeld non-linear QED 



D. K. Sinclair 



1. Introduction 

The n + 1 dimensional Born-Infeld (non-linear electrodynamics) action [jj], ||] is: 



S = b 2 J d n+l x 



1 - W-det I g^ v + -F^ v 



(1.1) 



This has seen a revival as a theory of strings and branes |(| [7p. Choosing « = 9, and 

dimensionally reducing this action from 9 + 1 to p + 1 dimensions describes a /?-brane. The 9 — /? 
additional components of are identified with the transverse components of the string/brane. 

Most of the serious work on these theories has dealt with their classical behaviour [§J. We 
are simulating these quantum theories on the lattice. We are starting with the simplest case where 
n = p = 3, the original Born-Infeld modification of electrodynamics, designed to make the self 
energy of a point charge finite. 

Section 2 reviews the classical Born-Infeld theory. In section 3 we indicate how this is ported 
to the lattice allowing Monte-Carlo simulations of the quantum theory. Section 4 details our simu- 
lations and preliminary results. A discussion of our results and conclusions are given in section 5. 



2. Classical Born-Infeld electrodynamics in Minkowski space-time 



This section summarises those results in references []9j, |1C|, |11|, [12p which are relevant for our 
investigations. 



Evaluating the determinant in equation LI the Lagrangian in 3 + 1 dimensions is 

^ = b 2 [l- 



1 



2r 



-b-\E-B) 2 }. 



(2.1) 



One can now define D and H by 



D 



H 



(EB)B 



dE 
dB 



l-b- 2 {E 2 -B 2 )-b-\E-B) 2 

B-b- 2 (E-B)E 
\-b- 2 {E 2 -B 2 )-b-\E-BY 



(2.2) 



Interaction with charged particles is implemented, as usual, by adding a term j^A 11 to the La- 
grangian. In terms of E, B, D and H, the equations of motion are the standard Maxwell equations. 
The non-linearity is hidden in equations |2.2[ 

For a static point charge p = e8 3 (r) the electric fields are 

e , 



D 
E 



Anr 2 
e f 



An 



r* + r* 



(2.3) 



where r = - r and ro = y Hence the D field for a static point charge is identical to the Maxwell 
solution, while the E field is screened at short distances. 
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3. Lattice Born-Infeld quantum-electrodynamics 

The Euclidean space action for Born-Infeld QED 

S = b 2 j d 4 x[^J 1 + b- 2 (E 2 + B 2 ) + b-\E • B) 2 - 1] 

is positive. Hence it can be simulated using Monte-Carlo methods. 
On the lattice we use the non-compact formulation: 

F^x+lfL + jV) = A v (x + fl)-A v (x)-A^(x + v)+A^(x) 



(3.1) 



(3.2) 



and average over the 16 choices of 6 plaquettes associated with each lattice site. We also define 
/3 = b 2 a 4 where a is the lattice spacing. Simulations are performed using the Metropolis Monte- 
Carlo method [13]. 

We measure the E and D fields due to a static point charge. This point charge e is introduced 
by including a Wilson Line (Polyakov Loop) W(x). 



W(x) = exp < ie^ 



A 4 (x,t) 



1 



Y AAy.t) 
N x N y N z y 4VJ ' ; 



(3.3) 



The second ('Jellium') term is needed, since a net charge would be inconsistent with periodic 
boundary conditions on A^. (E) and (D) in the presence of this charge are given by 



i(E) p (y-x) 
i(D) P (y-x) 



(E(y,t)W(x)) 

(W(x)) 
(D(y.t)W(x)) 

(W(x)) 



(3.4) 



Since W is complex, there is a sign problem, which causes (W(x)} to fall exponentially with N t . 
We use the method of Liischer and Weisz [[14]] (Parisi, Petronzio and Rapuano JT5|]) with thickness 
1 and 2 timeslices to overcome this exponential factor. 



4. Simulations and Results 



We have performed preliminary simulations of 500,000 10-hit Metropolis sweeps of the lattice 
at j8 = 100,1.0,0.01,0.0001 and 100,000 sweeps at j8 = 5,2,0.5,0.2,0.1, making measurements 
every 100 sweeps. We measured the E and D fields for on axis separations from the point charge. 
Figure [l] shows the expectation value of the Wilson lines (Polyakov loops) obtained from these 
simulations. Note that the value falls rapidly with increasing e. The rate of falloff also increases 
with increasing non-linearity (decreasing j8). Note also the small relative errors, even when the 
magnitude has fallen 8 orders of magnitude, which shows the effectiveness of the Liischer- Weisz 
method. 

Figure || shows the ratio of the E field in the direction of the separation from the charge to the 
D = E field for the free field (Maxwell) theory, at the minimum separation (Z = 0.5), in the limit of 
zero charge. At large j8, where the Born-Infeld theory asymptotes to the Maxwell theory, this ratio 
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Born — Infeld 8* Lattice 



Born— Infeld 8 Lattice 
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Figure 1: Wilson Lines as functions of charge e, for Figure 2: E /D ratio at minimum separation for e 
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as a function of j3 . 
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Figure 3: a) D fields at distance Z from the point charge <?, scaled by the free field (Maxwell) values, for 
j3 = 0.0001. b) E fields at distance Z from the point charge e, scaled by the free field (Maxwell) values, for 
/3 =0.0001. 



approaches 1. Classically, this ratio is 1 for all values of b. Quantum fluctuations cause this ratio to 
increase with increasing non-linearity (decreasing /3), approaching a value close to 3 for small j3. 

In figure |3|a, we plot the D fields scaled by their free field values at each separation, as a 
function of charge, for j8 = 0.0001, where the non-linearity is large, and the D field comes almost 
entirely from the (E -B)B term in its definition. The fact that this ratio is still 1 for all e at minimal 
separation is because D still obeys V - D = p, combined with cubic symmetry. Because we do not 
have full rotational symmetry on the lattice V • D = p is insufficient to make this ratio 1 at other 
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separations. The fact that this ratio is never more than 15% from 1 suggests that it would be 1 if we 
had rotational symmetry. However, we would expect rotational symmetry to be restored at large 
distances, which is why the ratio is closer to 1 for larger separations. Figure ||b is a similar graph 
for the E field. As well as showing the effects of quantum fluctuations as in figure ||, the E field is 
clearly screened at short distances as e increases, similar to what is seen in the classical theory. 

What is different from the classical theory is that classically the screening length continues to 
increase with decreasing b. The quantum theory approaches a limit as j3 — > 0. 



5. Discussion and conclusions 



We have succeeded in using lattice Monte-Carlo methods to extract non-perturbative physics 
from Born-Infeld electrodynamics, quantized using the Euclidean-time functional integral approach. 
The on-axis (quantum) electrostatic fields of a point charge are measured as functions of the charge 
e introduced as a Wilson Line. The approach of Liischer and Weisz, which reduces these mea- 
surements from an exponential- to a polynomial-time problem, was essential for extracting these 
quantities. 

In the classical field-theory E/D —> 1 as e —>0. For the quantum theory E /D increases from 1 
as the nonlinearity is increased indicating that the dielectric constant £ < 1. 

As \e\ is increased, the E field is screened at short distances. Screening increases with \e\ and 
with increasing nonlinearity (/3 ov b). The screening length ro appears to increase as y/\e\ as for 
the classical theory. D shows no such screening and appears to independent of the non-linearity, 
while D/e appears independent of e. 

Unlike the classical theory, where the screening length diverges as b — ► 0, the quantum theory 
approaches a fixed-point field theory as j8 = b 2 a 4 — ► 0. This conformal field theory has Euclidean 
Lagrangian Sf E = l\E.B\ and Hamiltonian Ji?=\DxB\ [|,[T|]. 



Normally, Born-Infeld QED is considered as an effective field theory with a momentum cutoff. 
However, as this cutoff — > oo it approaches the above fixed-point theory. If this fixed-point field 
theory is non-trivial, it would serve to define Born-Infeld QED without a cutoff. 

These first simulations were performed on 8 4 lattices. We are extending these simulations to 
larger lattices. We then plan to study those p-bmne theories obtained by dimensional reduction 
of n + 1 dimensional Born-Infeld theories to determine if the quantized theories continue to show 
string/brane dynamics. 

Acknowledgements 

Our simulations are currently running on the Rachael supercomputer at the Pittsburgh Super- 
computer Center. 

References 

[1] M. Born, On the Quantum Theory of Electromagnetic Field Proc. Roy. Soc. Lond. A 143 (1934) 410. 

[2] M. Born and L. Infeld, Foundations Of The New Field Theory, Proc. Roy. Soc. Lond. A 144 (1934) 
425. 



261/5 



Lattice simulations of Born-Infeld non-linear QED 



D. K. Sinclair 



[3] E. S. Fradkin and A. A. Tseytlin, Nonlinear Electrodynamics From Quantized Strings, Phys. Lett. B 
163(1985) 123. 

[4] R. G. Leigh, Dirac- Born-Infeld Action From Dirichlet Sigma Model, Mod. Phys. Lett. A 4 (1989) 
2767. 

[5] M. Aganagic, C. Popescu and J. H. Schwarz, Gauge-invariant and gauge-fixed D-brane actions, Nucl. 
Phys. B 495 (1997) 99 [arXiv:hep-th/96 12080]. 

[6] G. W. Gibbons, Born-Infeld particles and Dirichlet p-branes, Nucl. Phys. B 514 (1998) 603 
[arXiv:hep-th/9709027]. 

[7] A. A. Tseytlin, Born-Infeld action, supersymmetry and string theory, [arXiv:hep-th/9908105]. 

[8] C. G. Callan and J. M. Maldacena, Brane dynamics from the Born-Infeld action, Nucl. Phys. B 513 
(1998) 198 [arXiv:hep-th/9708147]. 

[9] I. Bialynicki-Birula, Nonlinear Electrodynamics: Variations On A Theme By Born And Infeld, In 
*Jancewicz, B. ( Ed.), Lukierski, J. ( Ed.): Quantum Theory Of Particles and Fields*, World 
Scientific(1983) 31-48. 

[10] D. Chruscinski and J. Kijowski, Generation of a dipole moment by external field in Born-Infeld 
non-linear electrodynamics, Compt. Rend. Acad. Sci. (Ser. II) 324 (1997) 435 
[arXiv:hep-th/9712101]. 

[11] D. Chruscinski, Point charge in the Born-Infeld electrodynamics, Phys. Lett. A 240 (1998) 8 
[arXiv:hep-th/97 12161]. 

[12] D. Chruscinski, Canonical formalism for the Born-Infeld particle, J. Phys. A 31 (1998) 5775 
[arXiv:hep-th/9803100]. 

[13] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller and E. Teller, Equation Of State 
Calculations By Fast Computing Machines, J. Chem. Phys. 21 (1953) 1087. 

[14] M. Luscher and P. Weisz, Locality and exponential error reduction in numerical lattice gauge theory, 
JHEP 0109 (2001) 010 [arXiv:hep-lat/0108014]. 

[15] G. Parisi, R. Petronzio and F. Rapuano, A Measurement Of The String Tension Near The Continuum 
Limit, Phys. Lett. B 128 (1983) 418. 

[16] D. Chruscinski, Strong field limit of the Born-Infeld p-form electrodynamics, Phys. Rev. D 62 (2000) 
105007 [arXiv:hep-th/0005215]. 



261/6 



